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ABSTRACT 

The Galaxy And Mass Assembly (GAMA) survey has been operating since February 
2008 on the 3.9-m Anglo- Australian Telescope using the AAOmega fibre-fed spec- 
trograph facility to acquire spectra with a resolution of R « 1300 for 120 862 SDSS 
selected galaxies. The target catalogue constitutes three contiguous equatorial regions 
centred at 9' 1 (G09), 12 h (G12) and 14.5' 1 (G15) each of 12 x 4 deg 2 to limiting fluxes 
of r pe t < 19.4, r pct < 19.8, and r pe t < 19.4 mag respectively (and additional limits 
at other wavelengths). Spectra and reliable rcdshifts have been acquired for over 98 
per cent of the galaxies within these limits. Here we present the survey footprint, 
progression, data reduction, redshifting, re-redshifting, an assessment of data quality 
after 3 years, additional image analysis products (including ugrizYJHK photome- 
try, Sersic profiles and photometric redshifts), observing mask, and construction of 
our core survey catalogue (GamaCore). From this we create three science ready cata- 
logues: GamaCoreDRl for public release, which includes data acquired during year 1 of 
operations within specified magnitude limits (February 2008 to April 2008); GamaCore- 
MainSurvey containing all data above our survey limits for use by the GAMA team 
© 2010 RAS and collaborators; and GamaCoreAtlasSv containing year 1, 2 and 3 data matched 

to Herschel- ATLAS Science Demonstration data. These catalogues along with the 
associated spectra, stamps, and profiles, can be accessed via the GAMA website: 
http : / / www . gama-survey . org/ 
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1 INTRODUCTION 

Large scale surveys are now a familiar part of the astronomy 
landscape and assist in facilitating a wide range of science 
programmes. Three of the most notable wide-area surveys 
in recent times, with a focus on galactic and galaxy evolu- 
tion, are the 2 Micron All Sky Survey (2MASS; Skrutskie et 
al. 2006), the 2 degree Field Galaxy Redshift Survey (2dF- 
GRS; Colless et al. 2001; 2003), and the Sloan Digital Sky 
Survey (SDSS; York et al. 2000). These surveys have each 
been responsible for a wide range of science advances at- 
tested by their publication and citation records (Trimble & 
Ceja 2010), ranging from: the identification of new stellar 
types (e.g., Kirkpatrick et al., 1999); tidal streams in the 
Galactic Halo (e.g., Belokurov et al. 2006); new populations 
of dwarf galaxies (e.g., Willman et al. 2005); galaxy popu- 
lation statistics (e.g., Bell et al. 2003; Baldry et al. 2006); 
the recent cosmic star-formation history (e.g., Heavens et 
al. 2004); group catalogues (e.g., Eke et al. 2004); merger 
rates (e.g., Bell et al., 2006); quantification of large scale 
structure (e.g., Percival et al. 2001); galaxy clustering (e.g., 
Norberg et al. 2001); and, in conjunction with CMB and 
SNIa searches, convergence towards the basic cosmological 
model now adopted as standard (e.g., Spergel et al. 2003; 
Cole et al. 2005). In addition to these mega-surveys there 
have been a series of smaller, more specialised, local sur- 
veys including the Millennium Galaxy Catalogue (MGC; 
Driver et al., 2005), the 6dF Galaxy Survey (6dFGS; Jones 
et al. 2004, 2009), the HI Parkes All Sky Survey (HIPASS; 
Meyer et al. 2004), and the Galaxy Evolution Explorer 
(GALEX; Martin et al. 2005) mission, which are each open- 
ing up new avenues of extra-galactic exploration (i.e., struc- 
tural properties, the near-IR domain, the 21cm domain, and 
the UV domain respectively). Together these surveys pro- 
vide an inhomogeneous nearby reference point for the very 
narrow high-z pencil beam surveys underway (i.e., DEEP2, 
VVDS, COSMOS, GEMS etc.), and from which compara- 
tive studies can be made to quantify the process of galaxy 
evolution (e.g., Cameron & Driver 2007). 

The Galaxy And Mass Assembly Survey (GAMA) has 
been established with two main sets of aims, which are sur- 
veyed in Driver et al. (2009). The first is to use the galaxy 
distribution to conduct a series of tests of the Cold Dark 
Matter (CDM) paradigm, and the second is to carry out de- 
tailed studies of the internal structure and evolution of the 
galaxies themselves. 

The CDM model is now the standard means by which 
data relevant to galaxy formation and evolution are inter- 
preted, and it has met with great success on 10-100 Mpc 
scales. The next challenge in validating this standard model 
is to move beyond small linear fluctuations, into the regime 
dominated by dark-matter haloes. The fragmentation of the 
dark matter into these roughly spherical virialized objects 
is robustly predicted both numerically and analytically over 
seven orders of magnitude in halo mass (e.g., Springel et al. 
2005). Massive haloes are readily identified as rich clusters 
of galaxies, but it remains a challenge to probe further down 
the mass function. For this purpose, one needs to identify 
low-mass groups of galaxies, requiring a survey that probes 
far down the galaxy luminosity function over a large rep- 
resentative volume. But having found low-mass haloes, the 
galaxy population within each halo depends critically on the 



interaction between the baryon processes (i.e., star forma- 
tion rate and feedback efficiency) and the total halo mass. 
In fact the ratio of stellar mass to halo mass is predicted 
(Bower et al. 2006; De Lucia et al. 2006) to be strongly 
dependent on halo mass, exhibiting a characteristic dip at 
Local Group masses. The need for feedback mechanisms to 
suppress star formation in both low mass haloes (via su- 
pernovae) and high mass haloes (via AGN) is now part of 
standard prescriptions in modeling galaxy formation. With 
GAMA we can connect these theoretical ingredients directly 
with observational measurements. 

However, the astrophysics of galaxy bias is not the only 
poorly understood area in the CDM model. Existing suc- 
cesses have been achieved at the price of introducing Dark 
Energy as the majority constituent of the Universe, and a 
key task for cosmology is to discriminate between various 
explanations for this phenomenon: a cosmological constant, 
time varying scalar field, or a deficiency in our gravity model. 
These aspects can be probed by GAMA in two distinct 
ways: either the form or evolution of the halo mass func- 
tion may diverge from standard predictions of gravitational 
collapse in the highly nonlinear regime, or information on 
non-standard models may be obtained from velocity fields on 
10-Mpc scales. The latter induce redshift-space anisotropies 
in the clustering pattern, which measure the growth rate 
of cosmic structure (e.g., Guzzo et al. 2008). Thus GAMA 
has the potential to illuminate both the astrophysical and 
fundamental aspects of the CDM model. 

Moving beyond the large-scale distribution, GAMA's 
main long-term legacy will be to create a uniform galaxy 
database, which builds on earlier local surveys in a compre- 
hensive manner to fainter flux levels, higher redshift, higher 
spatial resolution, and spanning UV to radio wavelengths. 
The need for a combined homogeneous, multiwavelength, 
and spatially resolved study can be highlighted by three top- 
ical issues: 

1. Galaxy Structure. Galaxies are typically comprised of 
bulge and/or disc components that exhibit distinct prop- 
erties (dynamics, ages, metallicities, profiles, dust and gas 
content), indicating potentially distinct evolutionary paths 
(e.g., cold smooth and hot lumpy accretion; cf. Driver et 
al. 2006 or Cook et al. 2010). This is corroborated by the ex- 
istence of the many SMBH-bulge relations (see, for example, 
Novak, Faber & Dekel 2006), which firmly couples spheroid- 
only evolution with AGN history (Hopkins et al. 2006). A 
comprehensive insight into galaxy formation and evolution 
therefore demands consideration of the structural compo- 
nents requiring high spatial resolution imaging on ~ 1 kpc 
scales or better (e.g., Allen et al. 2006; Gadotti 2009). 

2. Dust attenuation. A recent spate of papers (Shao et 
al. 2007; Choi et al. 2007; Driver et al. 2007, 2008; Masters 
et al. 2010) have highlighted the severe impact of dust atten- 
uation on the measurement of basic galaxy properties (e.g., 
fluxes and sizes). In particular, dust attenuation is highly 
dependent on wavelength, inclination, and galaxy type with 
the possibility of some further dependence on environment. 
Constructing detailed models for the attenuation of stellar 
light by dust in galaxies and subsequent re-emission (e.g., 
Popescu et al. 2000), is intractable without extensive wave- 
length coverage extending from the UV through to the far- 
IR. To survey the dust content for a significant sample of 
galaxies therefore demands a multi-wavelength dataset ex- 
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tending over a sufficiently large volume to span all environ- 
ments and galaxy types. The GAMA regions are or will be 
surveyed by the broader GALEX Medium Imaging Survey 
and Herschel-ATLAS (Eales et al. 2010) programmes, pro- 
viding UV to FIR coverage for a significant fraction of our 
survey area. 

3. The HI content. As star-formation is ultimately driven 
by a galaxy's HI content any model of galaxy formation and 
evolution must be consistent with the observed HI properties 
(see for example discussion in Hopkins, McClure-Griffiths & 
Gaensler 2008). Until recently probing HI beyond very low 
redshifts has been laborious if not impossible due to ground- 
based interference and/or sensitivity limitations (see for ex- 
ample Lah et al. 2009). The new generation of radio arrays 
and receivers are using radio-frequency interference mitiga- 
tion methods coupled with new technology receivers to open 
up the HI Universe at all redshifts (i.e., ASKAP, MeerKAT, 
LOFAR, and ultimately the SKA). This will enable coherent 
radio surveys that are well matched in terms of sensitivity 
and resolution to optical/near-IR data. Initial design study 
investment has been made in the DINGO project, which 
aims to conduct deep HI observations within a significant 
fraction of the GAMA regions using ASKAP (Johnston et 
al. 2007). 

The GAMA survey will eventually provide a wide- 
area highly complete spectroscopic survey of over 
400k galaxies with sub-arcsecond optical/near-IR imag- 
ing (SDSS/UKIDSS/VST/VISTA). Complementary multi- 
wavelength photometry from the UV (GALEX), mid- 
IR (WISE) and far-IR (Herschel) and radio wavelengths 
(ASKAP, GMRT) is being obtained by a number of inde- 
pendent public and private survey programmes. These ad- 
ditional data will ultimately be ingested into the GAMA 
database as they becomes available to the GAMA team. At 
the heart of the survey is the 3.9-m Anglo Australian Tele- 
scope (AAT), which is being used to provide the vital dis- 
tance information for all galaxies above well specified flux, 
size and isophotal detection limits. In addition the spec- 
tra from AAOmega for many of the sample will be of suffi- 
ciently high signal-to-noise and spectral resolution (~ 3 — 
6 A) to allow for the extraction of line diagnostic informa- 
tion leading to constraints on star-formation rates, velocity 
dispersions and other formation/evolutionary markers. The 
area and depth of GAMA compared to other notable surveys 
is detailed in Baldry et al. (2010; their Fig. 1). In general 
GAMA lies in the parameter space between that occupied 
by the very wide shallow surveys and the deep pencil beam 
surveys and is optimised to study structure on ~ 10 h^ 1 Mpc 
to 1 h^ 1 kpc scales, as well as sample the galaxy population 
from FUV to radio wavelengths. 

This paper describes the first three years of the GAMA 
AAOmega spectroscopic campaign, which has resulted in 
112k new redshifts (in addition to the 19k already known 
in these regions). In section 2 we describe the spectroscopic 
progress, data reduction, redshifting, re-redshifting, an as- 
sessment of the redshift accuracy and blunder rate, and an 
update to our initial visual classifications. In section 3 we 
describe additional image analysis resulting in ugrizY J H K 
matched aperture photometry, Sersic profiles, and photo- 
metric redshifts. In section 4 we describe the combination of 
the data presented in section 2 and 3 to form our core cat- 
alogue and investigate the completeness versus magnitude, 



colour, surface brightness, concentration, and close pairs. In 
section 5 we present the survey masks required for spatial 
clustering studies and in section 6 present our publicly avail- 
able science ready catalogue. These catalogues along with 
an MYSQL tool and other data inspection tools are now 
available at: http://www.gama-survey.org/ and we expect 
future releases of redshifts and other data products to occur 
on an approximately annual cycle. 

Please note all magnitudes used in this paper, unless 
otherwise specified, are r-band Petrosian (r pot ) from SDSS 
DR6, which have been extinction corrected and placed onto 
the true AB scale following the prescription described by the 
SDSS DR6 release (Adelman- McCarthy et al. 2008). 



2 THE GAMA AAT SPECTROSCOPIC 
SURVEY 

2.1 GAMA field selection, input catalogue and 
tiling algorithm 

The initial GAMA survey consists of three equatorial re- 
gions, each of 12 x 4 deg 2 (see Table [T] and Fig. [Q). The 
decision for this configuration was driven by three consider- 
ations: suitability for large-scale structure studies demand- 
ing contiguous regions of ~ 50 deg 2 to fully sample ~ 100 
co-moving h~ 1 Mpc structures at z ~ 0.2; observability de- 
manding a 5 hr Right Ascension baseline to fill a night's 
worth of observations over several lunations; and overlap 
with existing and planned surveys, in particular the SDSS 
(York et al. 2000), UKIDSS LAS (Lawrence et al. 2007), 
VST KIDS, VISTA VIKING, H-ATLAS (Eales et al. 2010), 
and ASKAP DINGO. Fig. [1] shows the overlap of some of 
these surveys. Note the H-ATLAS SGP survey region has 
changed since Fig. 5 of Driver et al. (2009) due to additional 
space-craft limitations introduced in-flight. The depth and 
area of the GAMA AAT spectroscopic survey were optimised 
following detailed simulations of the GAMA primary science 
goal of measuring the Halo Mass Function. This resulted in 
an initial survey area of 144 deg 2 to a depth of r pct < 19.4 
mag m the 9 h and 15 h regions and an increased depth of 
J" P ct < 19.8 mag in the 12 h region (see Table [TJ. In addition 
for year 2 and 3 observations additional K and z-band selec- 
tion was introduced such that the final main galaxy sample 
(Main Survey) can be defined (see Baldry et al., 2010) as 
follows: 

G09: r pct < 19.4 OR (if K ro n < 17.6 AND r modcl < 20.5) OR 
(zmodd < 18.2 AND r modc i < 20.5) mag 

G12: r p8 t < 19.8 OR (if K ron < 17.6 AND r modcl < 20.5) OR 
(Zmodoi < 18.2 AND r modc i < 20.5) mag 

G15: r pct < 19.4 OR (A' Kr on < 17.6 AND r modc i < 20.5) OR 
(Zmodei < 18.2 AND r modc i < 20.5) mag 

with all magnitudes expressed in AB. For the remainder 
of this paper we mainly focus, for clarity, on the r-band 
selected data and note that equivalent diagnostic plots to 
those shown later in this paper can easily be created for the 
z and K selections. 
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Table 1. Coordinates of the three GAMA equatorial fields. 



Field 


RA(J2000 deg.) 


Dcc(J2000 deg.) 


Area (deg 2 ) 


Depth (mag) 


G09 


129.0.. .141.0 


+3.0... - 1.0 


12 X 4 


r pot < 19.4 


G12 


174.0.. .186.0 


+2.0... - 2.0 


12 X 4 


r pot < 19.8 


G15 


211. 5. ..223. 5 


+2.0... - 2.0 


12 X 4 


r pet < 19.4 



-B!>- 




Figure 1. GAMA Phase I (black squares) in relation to other recent and planned surveys (see key). For a zoom in to the GAMA regions 
showing SDSS, UKIDSS and GALEX overlap please see Fig. 1 of the companion paper describing the photometry by Hill et al., (2010a). 
Also overlaid as grey dots are all known redshifts at z < 0.1 taken from NASA ExtraGalactic Database. 



2.2 Survey preparation 

The input catalogue for the GAMA spectroscopic survey 
was constructed from SDSS DR6 (Adelman-McCarthy et 
al. 2008) and our own reanalysis of UKIDSS LAS DR4 
(Lawrence et al. 2007) to assist in star-galaxy separation; see 
Hill et al. (2010a) for details. The preparation of the input 
catalogue, including extensive visual checks and the revised 
star-galaxy separation algorithm, is described and assessed 
in detail by Baldry et al. (2010). The survey is being con- 
ducted using the AAT's AAOmega spectrograph system (an 



upgrade of the original 2dF spectrographs: Lewis et al. 2002; 
Sharp et al. 2006). 

In year 1 we implemented a uniform grid tiling algo- 
rithm (see Fig. [2] upper panels) , which created a significant 
imprint of the tile positions on the spatial completeness dis- 
tribution. Subsequently, for years 2 and 3 we implemented 
a heuristic "greedy" tiling strategy, which was designed to 
maximise the spatial completeness across the survey regions 
within 0.14 deg smoothed regions. Full details of the GAMA 
science requirements and the tiling strategy devised to meet 
these are laid out by Robotham et al. (2010). The efficiency 
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of this strategy is discussed further in section [S] The final 
location of all tiles is shown in Fig. [2] and the location of 
objects for which redshifts were not secured or not observed 
are shown in the lower panels of Fig. [3] as black dots or red 
crosses respectively. 



2.3 Observations and data reduction 

All GAMA 2dF pointings (tiles) were observed during dark 
or grey time with exposure times mostly ranging from 3000 
to 5000 s (in 3 to 5 exposures) depending on seeing and 
sky brightness. Observations were generally conducted at 
an hour angle of less than 2 hr (the median zenith dis- 
tance of the observations is 35°) and with the Atmospheric 
Dispersion Corrector engaged. We used the 580V and 385R 
gratings with central wavelengths of 4800 A and 7250 A in 
the blue and red arms, respectively, separated by a 5700 A 
dichroic. This set-up yielded a continuous wavelength cov- 
erage of 3720A-8850A at a resolution of ~ 3.5 A (in the 
blue channel) and « 5.5 A (in the red channel). Every block 
of science exposures was accompanied by a flat-field and 
an arc-lamp exposure, while a master bias frame was con- 
structed once per observing run. In each observation fibres 
were allocated to between 318 and 366 galaxy targets (de- 
pending on the number of broken fibres), to at least 20 
blank sky positions, to 3-5 SDSS spectroscopic standards, 
and to 6 guide stars selected from the SDSS in the range 
14.35 < r-p t < 14.5 mag (verified via cross-matching with 
the USNO-B catalogue to ensure proper motions were below 
15 mas/yr). Galaxy targets were prioritised as indicated in 
Table 1 of Robotham et al. (2010). The location of all 392 
tile centres are shown in Fig. [2] The change in tiling strat- 
egy from a fixed grid system in year 1 to the greedy tiling 
strategy in years 2 and 3 is evident with the latter leading to 
an extremely spatially uniform survey (see section [5]). The 
overall progress of the survey in terms of objects per night 
and the cumulative total numbers are shown in Fig. [3] illus- 
trating that typically between 1500 and 2500 redshifts were 
obtained per night over the three year campaign. 

The data were reduced at the telescope in real time 
using the (formeiQ) Anglo-Australian Observatory's 2dfdr 
software (Croom, Saunders & Heald 2004) developed con- 
tinuously since the advent of 2dF, and recently optimised 
for AAOmega. The software is described in a number of 
AAO documentfl Briefly, it performs automated tramline 
detection, sky subtraction, wavelength calibration, stacking 
and splicing. Following the standard data reduction we at- 
tempted to improve the sky subtraction by applying a prin- 
cipal component analysis (PCA) technique similar to that 
described by Wild & Hewett (2005) and described in detail 
by Sharp & Parkinson (2010). We note that ~ 5 per cent of 
all spectra are affected by fringing (caused by small air gaps 
between the adhesive that joins a fibre with its prism), and 
we are currently studying algorithms that might remove or 
at least mitigate this effect. 



1 Note that from 1st July 2010 the Anglo Australian Observatory 
(AAO) has been renamed the Australian Astronomical Observa- 
tory (AAO), see Watson & Colless (2010). 
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Yearl (Feb-May) Year 2 (Feb-May) Year 3 (Feb-May) 

Figure 4. Progression of the GAMA survey in terms of spectra 
and redshift acquisition for years 1, 2 and 3 as indicated. The 
cyan line shows the number of spectra obtained per night within 
our main survey limits, mauve shows the number below our flux 
limits (i.e., secondary targets and fillers). The solid line shows the 
cumulative distribution of redshifts within our survey limits and 
the dashed line shows the cumulative distribution of all redshifts. 
The cumulative distributions include pre-existing redshifts and 
are calculated monthly rather than nightly. 



2.4 Redshift ing 

The fully reduced and PCA-sky-subtracted spectra were ini- 
tially redshifted by the observers at the telescope using the 
code runz, which was originally developed by Will Suther- 
land for the 2dFGRS (now maintained by Scott Croom). 
runz attempts to determine a spectrum's redshift (i) by 
cross-correlating it with a range of templates, including star- 
forming, E+A and quiescent galaxies; A, O and M stars; as 
well as QSO templates; (ii) by fitting Gaussians to emission 
lines and searching for multi-line matches. These estimates 
are quasi-independent because the strongest emission lines 
are clipped from the templates before the cross-correlation 
is performed, runz then proceeds by presenting its operator 
with a plot of the spectrum marking the positions of com- 
mon nebular emission and stellar absorption lines at the best 
automatic redshift. This redshift is then checked visually by 
the operator who, if it is deemed incorrect, may use a num- 
ber of methods to try to find the correct one. The process is 
concluded by the operator assigning a (subjective) quality 
(Q) to the finally chosen redshift: 

Q = 4: The redshift is certainly correct. 

Q = 3: The redshift is probably correct. 

Q = 2: The redshift may be correct. Must be checked before 

being included in scientific analysis. 

Q = 1: No redshift could be found. 

Q = 0: Complete data reduction failure. 

With the above definitions it is understood that by assigning 
Q > 3 a redshift is approved as suitable for inclusion in 
scientific analysis. Note that Q refers to the (subjective) 
quality of the redshift, not of the spectrum. 



2.5 Re-redshifting 

The outcome of the redshifting process described above will 
not be 100 per cent accurate. It is inevitable that some frac- 
tion of the Q > 3 redshifts will be incorrect. Furthermore, 



http://www.aao.gov.au/AAO/2df/aaomcga/aaomega_software.html the quality assigned to a redshift is somewhat subjective and 
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Figure 2. Location of the tiles in year 1 (top), year 2 (centre), and year 3 (bottom) and for G09 (left), G12 (centre) and G15 (right). 



will depend on the experience of the redshifter. In an effort 
to weed out mistakes, to quantify the probability of a red- 
shift being correct, and thereby to homogenise the quality 
scale of our redshifts, a significant fraction of our sample 
has been independently re-redshifted. This process and its 
results will be described in detail by Liske et al. (2010, in 
prep.). Briefly, the spectra of all Q = 2 and 3 redshifts, 
of all Q — 4 redshifts with discrepant photo-z's, and of an 
additional random Q = 4 sample have been independently 
re-redshifted, where those conducting the re-redshifting had 
no knowledge of the originally assigned redshift or Q. The 
Q = 2 sample was re-examined twice. Overall approximately 
one third of the entire GAM A sample was re-evaluated. The 
results of the blind re-redshifting process were used to esti- 
mate the probability, for each redshifter, that she/he finds 
the correct redshift as a function of Q (for Q > 2), or that 
she/he has correctly assigned Q = 1. Given these probabili- 
ties, and given the set of redshift "opinions" for a spectrum, 
we have calculated for each redshift found for this spectrum 
the probability, p z , that it is correct. This allowed us to se- 
lect the "best" redshift in cases where more than one redshift 
had been found for a given spectrum. It also allowed us to 
construct a "normalised" quality scale: 

nQ = 4 if p z > 0.95 

nQ = 3 if 0.9 < p z < 0.95 

nQ = 2 if p z < 0.9 

nQ = 1 if it is not possible to measure a redshift from this 
spectrum. 

Unlike Q, whose precise, quantitative meaning depends on 
the redshifter who assigned it, the meaning of nQ is homo- 
geneous across the entire redshift sample. 

In Fig. [5] we show one example each of a spectrum with 
nQ — 4, nQ = 3 and nQ = 2. 

2.6 Final redshift sample 

In the three years of observation completed so far we have 
observed 392 tiles, resulting in 135 902 spectra in total 



(including standard stars), and 134 390 spectra of 120 862 
unique galaxy targets. Of these, 114 043 have a reliable red- 
shift with nQ > 3, implying a mean overall redshift com- 
pleteness of 94.4 per cent. Restricting the sample to the 
r-limited Main Survey targets (i.e., ignoring the K and z 
selection and fainter fillers), we find that the completeness 
is > 98 per cent in all three GAMA regions, leaving little 
room for any severe spectroscopic bias. Fig.[6]shows the evo- 
lution of the survey completeness for the main r-band lim- 
ited sample versus apparent magnitude across (left-to-right) 
the three GAMA regions. 

2.7 Redshift accuracy and reliability 

Quantifying the redshift accuracy and blunder rate is cru- 
cial for most science applications and can be approached 
in a number of ways. Here we compare redshifts obtained 
for systems via repeat observations within GAMA (Intra- 
GAMA comparison) and also repeat observations of objects 
surveyed by earlier studies (Inter-survey comparison). 

2.7.1 Intra-GAMA comparison 

Our sample includes 974 objects that were observed more 
than once and for which we have more than one nQ > 3 
redshift from independent spectra. The distribution of pair- 
wise velocity differences, Av, of this sample is shown as the 
shaded histogram in the left panel of Fig. [7] This distribution 
is clearly not Gaussian but roughly Lorentzian (blue line), 
although with a narrower core, and there are a number of 
outliers (see below). Nevertheless, if we clip this distribution 
at ±500 km s _1 we find a 68-percentile range of 185 km s _ , 
indicating a redshift error o~ v = 65 km s~ . However, this 
value is likely to depend on nQ. Indeed, if we restrict the 
sample to pairs where both redshifts have nQ — 4 (red his- 
togram), we find cr„,4 = 60 km s _ . Our sample of pairs 
where both redshifts have nQ = 3 is small (22 pairs) but 
this yields a Vt 3 = 101 km s _ . However, given cr„,4 we can 
also use our larger sample of pairs where one redshift has 
nQ — 3 and the other has nQ = 4 (green histogram) to 
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Figure 3. (top) The distribution of pre-existing redshifts (SDSS, yellow; 2dFGRS, red; MGC, green; other, blue) within the three 
regions (as indicated), (middle) the distributions of redshifts acquired during the GAMA Phase I campaign, and (lower) the distribution 
of targets for which redshifts were not obtained (black dots) or were not targeted (red crosses). 



obtain an independent estimate of a Vt 3 = 97 km s , which 
is in reasonable agreement. 

Defining discrepant redshift pairs as those with \Av\ > 
500 km s _1 , and assuming that only one, but not both of 
the redshifts of such pairs is wrong, we find that 3.6 per 
cent of redshifts with nQ = 4 are in fact wrong. This is in 
reasonable agreement with the fact that nQ = 4 redshifts 
are defined as those with p z > 0.95. However, for nQ — 3 
we find a blunder rate of 15.1 per cent, which is somewhat 
higher than expected based on the fact that nQ = 3 is de- 
fined as p z > 0.9. We surmise that this is likely to be the 
result of a selection effect: many of the objects in this sample 
are likely to have been re-observed by GAMA because the 



initial redshift of the first spectrum was only of a low quality 
(i.e. Q = 2). However, subsequent re-redshifting of the ini- 
tial spectrum (after the re-observation) may have produced 
confirmation of the initial redshift, which will have bumped 
the redshift quality to nQ — 3. Hence this sample is likely 
to include many spectra that are of worse quality, and or 
spectra that are harder to redshift than the average nQ = 3 
spectrum, which will produce a higher blunder rate for this 
sample than the average blunder rate for the whole nQ = 3 
sample. Indeed, the median S/N of this sample is 20 per cent 
lower than for the full GAMA sample. Note that this may 
also have an effect on the redshift accuracies determined 
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Figure 5. Examples of spectra with redshift quality nQ = 4 (top panel), 3 (middle) and 2 (bottom). We show the spectrum (black), 
the 1(T error (green) and the mean sky spectrum (blue, scaled arbitrarily w.r.t. the spectrum). The vertical dashed red lines mark the 
positions of common nebular emission and stellar absorption lines at the redshift of the galaxy. The spectra were smoothed with a boxcar 
of width 5 pixels. 
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Figure 6. Evolution the redshift completeness (Q > 3) of the GAMA survey (main r-band selection only) over three years of observations 
showing the progressive build-up towards uniform high completeness. The horizontal line denotes a uniform 95 per cent completeness. 



above. This will be investigated in more detail by Liske et 
al. (2010, in prep). 



2.7.2 Inter-survey comparisons 

Our sample includes 2522 unique GAMA spectra (with 
nQ > 3 redshifts) for objects that had previously been ob- 
served by other surveys (see section [2~8)l . for a total of 2671 
GAMA-non-GAMA pairs. The distribution of the velocity 
differences of this sample is shown as the shaded histogram 
in the right panel of Fig. Approximately 81 per cent of 
the non-GAMA spectra in this sample are from the 2dFGRS 
and the MGC, which were both obtained with 2dF, using the 
same set-up and procedures. For the 2dFGRS, Colless et al. 
(2001) quote an average redshift uncertainty of 85 km s~ . 
Using this value together with the observed 68-percentile 
ranges of the velocity differences of the GAMA(nQ = 4, 3)- 
non-GAMA(nQ > 3) sample (red and green histograms, re- 
spectively) we find a Vl 4 = 51 km s~ and cr„,3 = 88 km s - . 
These values are somewhat lower than those derived in the 
previous section. The most likely explanation for this is that 
the inter-survey sample considered here is brighter than the 
intra-GAMA sample of the previous section (because of the 
spectroscopic limits of the 2dFGRS and MGC). Indeed, the 
median S/N of the GAMA spectra in the inter-survey sam- 
ple is a factor of 1.7 higher than that of the spectra in the 
intra-GAMA sample. 

As above, we can also attempt to estimate the GAMA 
blunder rate from the inter-survey sample. From the 
GAMA(nQ = 4)-non-GAMA(nQ > 3) sample we find a 
GAMA nQ = 4 blunder rate of 5.0 per cent if we assume 
that all redshift discrepancies are due to GAMA mistakes. 



However, this is clearly not the case since the blunder rate 
improves to 3.0 per cent if we restrict the sample to pairs 
with nQ > 4 non-GAMA redshifts. Similarly, the GAMA 
nQ = 3 blunder rate comes out at 10.1 or 4.2 per cent de- 
pending on whether one includes the pairs with non-GAMA 
nQ = 3 redshifts or not, considerably lower than the corre- 
sponding value derived in the previous section. 

In summary, it seems likely that neither the intra- 
GAMA nor the inter-survey sample are fully representative 
of the complete GAMA sample. The former is on average 
of lower quality than the full sample while the latter is of 
higher quality. Hence we must conclude that the redshift ac- 
curacies and blunder rates determined from these samples 
are not representative either, but are expected to span the 
true values. A full analysis of these issues will be provided 
by Liske et al. (2010) following re-redshifting of the recently 
acquired year 3 data. 



2.8 Merging GAMA with data from earlier 
redshift surveys 

The GAMA survey builds upon regions of sky already sam- 
pled by a number of surveys, most notably the SDSS and 
2dFGRS but also several others (as indicated in Table[2}. As 
the GAMA input catalogue (Baldry et al. 2010) has taken 
the pre-existing redshifts into account, the GAMA data by 
themselves constitute a highly biased sample missing 80- 
90 per cent of bright sources with r pct < 17.77 mag and 
fainter objects previously selected by AGN or LRG surveys. 
It is therefore important for almost any scientific applica- 
tion, outside of analysing AAOmega performance, to pro- 
duce extended catalogues that include both the GAMA and 
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Figure 7. Left: The shaded histogram shows the distribution of differences between the redshifts measured from independent GAMA 
spectra of the same objects, where all redshifts have nQ > 3 (868 pairs from 1718 unique spectra of 856 unique objects). The red 
histogram shows the same for pairs where both redshifts have nQ = 4 (617 pairs from 1216 unique spectra of 605 unique objects). The 
green histogram shows the same for pairs where one redshift has nQ = 3 and the other nQ = 4 (229 pairs from 458 unique spectra 
of 229 unique objects). The blue line shows a Lorentzian with 7 = 50 km s _1 for comparison. Right: The shaded histogram shows the 
distribution of differences between the redshifts measured from independent GAMA and non-GAMA spectra of the same objects, where 
all redshifts have nQ > 3 (2533 pairs from 4892 unique spectra of 2359 unique objects). The red histogram shows the same for pairs 
where the GAMA redshift has nQ = 4 and the non-GAMA redshift nQ > 3 (2385 pairs from 4618 unique spectra of 2233 unique objects). 
The green histogram shows the same for pairs where the GAMA redshift has nQ = 3 (148 pairs from 290 unique spectra of 142 unique 
objects). The blue line shows a Lorentzian with 7 = 70 km s _1 for comparison. 



pre-GAMA data. Fig.[8]shows the n(z) distributions of red- 
shifts in the three GAMA blocks to r pct < 19.4 mag in G09 
and G15 and to r pct < 19.8 mag in G12 colour coded to 
acknowledge the survey from which they originate. Table [2] 
shows the contribution to the combined redshift catalogue 
from the various surveys. Note the numbers shown in Ta- 
ble [2] may disagree with those shown in Baldry et al. (2010) 
as some repeat observations of previous targets were made. 



2.9 Update to visual inspection of the input 
catalogue 

The GAMA galaxy target catalogue has been constructed 
in an automated fashion from the SDSS DR6 (see Baldry et 
al. 2010) and a number of manual checks of the data based on 
flux and size ratios and various SDSS flags were made. This 
resulted in 552 potential targets being expunged from the 
survey prior to Year 3 (Baldry et al. 2010). Expunged targets 
were given VIS_CLASS values of 2 (no evidence of galaxy light) 
or 3 (not the main part of a galaxy). In reality, objects with 
VIS_CLASS=3 were targeted but at a lower priority (below the 
main survey but above any filler targets) with most receiving 
redshifts by the end of Year 3. About 6 percent of these were, 
by retrospective visual inspection and by the difference in 
redshift, clearly not part of the galaxy to which they were 
assigned. These 16 objects were added back into the Main 
Survey (vis_CLASS=l). 

After the redshifts were all assigned to the survey ob- 
jects, we further inspected two distinct categories: very 
bright objects (r pe t < 17.5 mag, 9453 objects); and tar- 
gets for which no redshift was recovered (2083 objects). The 
original visual classification was made using the SDSS jpg 
image tools (Lupton et al. 2004; Neito-Santisteban, Szalay 
& Gray 2004). Here we also created postage stamp images 



by combining the it, r, K images, and the resulting 11536 im- 
ages were visually inspected by SPD. It became clear that 
many of the missed apparently bright galaxies were in fact 
probably not galaxies or at least much fainter; and some of 
the other apparent targets were also probably much fainter. 
A VIS_CLASS value of 4 was introduced meaning "compro- 
mised photometry (selection mag has serious error)". From 
the inspection, 50 objects were classified as VIS_CLASS=4 
and 40 objects as VIS_CLASS=2. Independent confirmation 
was made by 1KB using the SDSS jpg image tools. In sum- 
mary, a total of 626 potential targets identified by the au- 
tomatic prescription described in Baldry et al. (2010) were 
expunged from the survey because of the visual classifica- 
tion (2 < VIS_CLASS < 4). These objects are identified in the 
GamaTiling catalogue available from the data release web- 
site. 

The final input catalogue (GamaTiling) therefore con- 
stitutes 30331, 50924, and 33261 objects above the Main 
r-band Survey limits of r pct < 19.4, r pct < 19.8, and 
r pct < 19.4 mag in G09, G12 and G15 respectively. Fig. [g] 
shows the normalised galaxy number-counts in these three 
fields indicating a significant variation between the three 
fields and in particular a significant underdensity in G09 to 
r pc t = 19.0 mag. This will be explored further in section [4731 



3 ADDITIONAL IMAGE ANALYSIS 
PRODUCTS 

At this point we have two distinct catalogues: the input 
catalogue (GamaTiling) used by the tiling algorithm, which 
defines all viable target galaxies within the GAMA regions; 
and the combined redshift catalogue (GamaRedshifts), which 
consists of the pre-existing redshifts and those acquired by 
the GAMA Team and described in detail in the previous 
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Figure 8. The n(z) distributions of the GAMA regions including new GAMA rcdshifts (shown in yellow) alongside pre-existing rcdshifts 
already in the public domain as indicated. 



Table 2. Contribution from various surveys to the final GAMA database. 



Survey 


G09 


G12 


G15 


Reference or 


source (Z.SOURCE) 


r < 19.4 


r < 19.8 


r < 19.4 


Acknowledgment 


SDSS DR7 


3190 


4758 


5092 


Abazajian ct al. (2009) 


2dFGRS 





2107 


1196 


Colless et al. (2001) 


MGC 





612 


497 


Driver et al. (2005) 


2SLAQ-LRG 


2 


49 


13 


Canon ct al. (2006) 


GAMAz 


26783 


42210 


25858 


This paper 


6dFGS 


7 


14 


10 


Jones ct al. (2009) 


UZC 


3 


2 


1 


Falco et al. (1999) 


2QZ 





31 


5 


Croom ct al. (2004) 


2SLAQ-QSO 


1 


1 


1 


Croom ct al. (2009) 


NED 





2 


3 


The NASA Extragalactic Database 


z not known 


347 


1145 


591 




Total 


30333 


50931 


33267 




Completeness (%) 


98.9% 


97.8% 


98.0 % 





section. Our core catalogues consist of the combination of 
these two catalogues with three other catalogues containing 
additional measurements based on the ugrizY J H K imaging 
data from SDSS DR6 and UKIDSS LAS archives. These ad- 
ditional catalogues are: GamaPhotometry, GamaSersic, and 
GamaPhotoz and are briefly described in the sections below 
(for full details see Hill et al, 2010a; Kelvin et al. in prep, 
and Parkinson et al. in prep) . These five catalogues are then 
combined and trimmed via direct name matching to produce 
the GAMA Core dataset (GamaCore), which represents the 
combined data from these five catalogues. 



3.1 ugrizYJHK matched aperture photometry 
(GamaPhotometry) 

A key science objective of GAMA is to provide cross- 
wavelength data from the far ultraviolet to radio wave- 
lengths. The first step in this process is to bring together 
the available optical and near-IR data from the SDSS and 
UKIDSS archives. Initially the GAMA team explored sim- 
ple table matching, however this highlighted a number of 
issues, most notably: differences in deblending outcomes, 
aperture sizes, and seeing between the SDSS and UKIDSS 
datasets. Furthermore as the UKIDSS YJHK data frames 
are obtained and processed independently in each band by 
The Cambridge Astronomical Survey Unit group there is 
the potential for inconsistent deblending outcomes, incon- 
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Figure 9. Normalised number-counts as indicated for the three 
regions and the average. Error bars are based on Poisson statistics 
suggesting an additional source of error, most likely cosmic (sam- 
ple) variance between the three GAMA fields. The data becomes 
consistent only at r pc t > 19 mag. 



sistent aperture sizes and seeing offsets between the YJHK 
bandpasses. Finally both SDSS and UKIDSS LAS measure 
their empirical magnitudes (Petrosian and Kron) using cir- 
cular apertures (see Kron 1980 and Petrosian 1976). For a 
partially resolved edge-on system this can result in either 
underestimating flux or adding unnecessary sky signal in- 
creasing the noise of the flux measurement and compound- 
ing the deblending issue. To overcome these problems the 
GAMA team elected to re-process all available data to pro- 
vide matched aperture photometry from it to K using ellip- 
tical Kron and Petrosian apertures defined using the r-band 
data (our primary selection waveband). This process and 
comparisons to the original data are outlined in detail in 
Hill et al. (2010a) and described here in brief as follows: 

All available data are downloaded and scaled to a uni- 
form zeropoint in the strict AB magnitude system (i.e., after 
first correcting for the known SDSS to AB offset in it and z 
bands and converting the UKIDSS LAS from Vega to AB, 
see Hill et al. 2010b for conversions). Each individual data 
frame is convolved using a Gaussian PSF to yield consistent 
2.0" FWHM PSF measurements for the intermediate flux 
stars. The data frames (over 12,000) are then stitched into 
single images using the SWARP software developed by the 
TERAPIX group (Bertin et al. 2002) resulting in twenty- 
seven 20 Gb images each of 12 x 4 deg 2 . at 0.4" x 0.4" pixel 
sampling. The SWARP process removes the background us- 
ing the method described in SExtractor (Bertin & Arnouts 
1996) using a 256 x 256 pixel mesh. SExtractor is then ap- 
plied to the r-band data to produce a master catalogue and 
rerun in each of the other 8 bandpasses in dual object mode 
to ensure consistent r-defined Kron apertures (2.5Rk, see 
Kron 1980) from u to K. As the data have been convolved 
to the same seeing no further correction to the colours are re- 



Figure 10. (Upper) all main survey galaxies are plotted accord- 
ing to their colour relative to J-band at their observed wave- 
length. The data points combine to make a global spectral en- 
ergy distribute for the galaxy population at large. Outlier points 
are typically caused by an erroneous data point in any particu- 
lar filter, (center) the same plot but using the revised photome- 
try. The lower panel shows the 5<x-clipped standard deviation for 
these two distributions indicating that in the uYJHK bands the 
colour distribution is quantifiably narrower reflecting an improve- 
ment in the photometry, providing a cleaner dataset for detailed 
SED modeling. Only data with photometry in all bands and with 
secure redshifts less than 0.5 are shown for clarity. 



quired. For full details see Hill et al., (2010a). Fig. [T0l shows 
the it to K data based on archival data (upper) and our 
reanalysis of the archival data (centre). The data points are 
plotted as a colour offset with respect to z-band at their rest 
wavelengths, hence the data stretch to lower wavelengths 
within each band with redshift. As the k and e corrections 
are relatively small in z-band the data dovetails quite nicely 
from one filter to the next. The lower panel shows the stan- 
dard deviation within each log wavelength interval indicat- 
ing the colour range over which the data are spread. In the 
U and YJHK bands the colour distribution is noticeably 
narrower for the revised photometry (blue line) over the 
archival photometry (red line) . This indicates a quantifiable 
improvement in the photometry in these bands over the orig- 
inal archival data with the majority of the gain occurring 
for objects in more crowded regions. GamaPhotometry con- 
tains r-band detections with matched aperture photometry 
for ~ 1.9 million objects, 1 million of which are matched to 
SDSS DR6 objects. The photometric pipeline is described in 
full in Hill et al. (2010a) and all SWARP images are made 
publicly available for downloading at the GAMA website. 

3.2 Sersic profiles (GamaSersic) 

By design both Petrosian and Kron magnitude systems re- 
cover only a proportion of a galaxy's flux. For the two most 
commonly discussed galaxy profiles, the exponential and the 
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de Vaucouleurs profiles, traditionally used to describe galaxy 
discs or spheroids, the Kron measurement recovers 96 and 
90 per cent while the SDSS implementation of the Petrosian 
profile recovers 98 per and 83 per cent respectively (see 
Graham & Driver 2005). These two profiles are two spe- 
cific cases of the more general Sersic profile introduced by 
Sersic (1963, 1968) and more recently reviewed by Graham 
& Driver (2005). The Sersic profile (see Eqn. [1} is a use- 
ful general description of a galaxy's overall light profile and 
has also been used to profile the dark matter distribution in 
numerical simulations (see Merritt et al. 2006): 

/(r) = / exp(-(r/a) 1/ ' 1 ). (1) 

The Sersic model has three primary parameters (see Eqn.Q]), 
which are the central intensity (7 ); the scale- length (a); the 
Sersic index (n); and two additional parameters for defining 
the ellipticity (e) and the position angle (8). Note this ex- 
pression can also be recast in terms of the effective surface 
brightness and half-light radius (see Graham & Driver 2005 
for a full description of the Sersic profile). The Sersic index 
might typically range from n — 0.5 for diffuse systems to 
n — 15 for concentrated systems. If a galaxy's Sersic pro- 
file is known it is possible to integrate this profile to obtain 
a total magnitude measurement, and this mechanism was 
used by the SDSS to provide model magnitudes by force 
fitting either an n = 1 (exponential) or an n = 4 (de Vau- 
couleurs) profile to all objects. Here we attempt to take the 
next step, which is to derive the Sersic profile with n as a 
free parameter, in order to provide total magnitudes for all 
targets. 

As described in full in Kelvin et al. (2010, in prep.) 
we use the GALFIT3 package (Peng et al. 2010) to fit all 
galaxies to SDSS DR6 r pct < 22.0 mag. In brief the pro- 
cess involves the construction of comparable SWARP'ed im- 
ages as in section [37X1 but without Gaussian PSF convolution 
(raw SWARPS). Instead, the PSF at the location of every 
galaxy is derived using 20 intermediate brightness stars from 
around the galaxy as identified by SExtractor and modeled 
by the code PSFex (Bertin priv. comm.). Using the model 
PSF the 2D Sersic profile is then derived via GALFIT3, 
which convolves a theoretical profile with the PSF and min- 
imises the five free parameters. The output Sersic total mag- 
nitude is an integral to infinite radius, which is unrealistic. 
However relatively little is known as to how the light profile 
of galaxies truncate at very faint isophotes with all variants 
seen (Pohlen & Trujillo 2006). 

In order to calculate an appropriate Sersic magnitude 
it is prudent to adopt a truncation radius out to which the 
Sersic profile is integrated. For SDSS model magnitudes ex- 
ponential (n = 1) profiles were smoothly truncated beyond 
3i? e and smoothly truncated beyond 7R e for de Vaucouleurs 
(n = 4) profiles. Here we adopt an abrupt truncation ra- 
dius of 10Ji e ; for the majority of our data, this equates to 
an isophotal detection limit in the r-band of ~ 30 mag/sq 
arcsec — the limit to which galaxy profiles have been ex- 
plored. We note for almost all plausible values of n, the 
choice of a truncation radius of 7R e or 10i? e introduces 
a systematic uncertainty that is comparable or less than 
the photometric error. Fig. [TTJ shows a comparison of the 
three principle photometric methods, with r pct — rKron versus 
log 10 (n) (upper), r pct - r S ersici0R o versus log 10 (n) (middle), 
and r K ron - rs cr sicioR c versus log 10 (n) (lower). Note that it 



is more logical to compare to log 10 (n) than n as it appears 
in the exponent of the Sersic profile definition (see review of 
the Sersic profile in Graham & Driver 2005 for more infor- 
mation). The top panel indicates that our Kron photometry 
(defined in section [3~Tj) . reproduces (with scatter) the origi- 
nal SDSS photometry. As described in the previous section 
the motivation for rederiving the photometry is to provide 
elliptical matched apertures in ugrizY JHK . The middle 
and lower plots shows a significant bias between these tradi- 
tional methods (Petrosian and Kron) and the new SersiclO 
R e magnitude. The vertical lines indicate the location of 
the canonical exponential and de Vaucouleurs profiles. For 
normal discs and low-n systems the difference is negligible. 
However for n > 4 systems the difference becomes signifi- 
cant. This is particularly crucial as the high-n systems are 
typically the most luminous and most massive galaxy sys- 
tems (Driver et al. 2006). Hence a bias against this popu- 
lation can lead to serious underestimation of the integrated 
properties of galaxies, e.g., integrated stellar mass or lumi- 
nosity density (see Fig. 21 of Hill et al. 2010a, where a 20 
per cent increase in the luminosity density is seen when mov- 
ing from Petrosian to Sersic magnitudes). The Sersic 10R e 
magnitudes are therefore doing precisely what is expected, 
which is to recover the flux lost via the Petrosian or Kron 
photometric systems. 

GamaSersic contains the five parameter Sersic output 
for all 1.2 million objects with SDSS DR6 r pct < 22 and 
full details of the process are described in Kelvin et al. (in 
prep.). At the present time we only include in our core 
catalogues the individual profiles for each object and the 
r_SERS_MAG_10RE -parameter (see Table |A2j) while final 
checks are being conducted. We release the total correction 
to allow for an early indication of the Kron- Total mag offset 
but caution that these values may change. Similarly we in- 
clude for individual objects all profiles with the full Sersic in- 
formation specified but again caution that these may change 
slightly over the next few months as further checks continue. 
If colour gradients are small then the r-band Kron to To- 
tal magnitude offset given by r_SERS_MAG_10RE can be 
applied to all filters. 

3.3 Photometric redshifts (GamaPhotoz) 

The GAMA spectroscopic survey is an ideal input dataset 
for creating a robust photometric redshift code that can then 
be applied across the full SDSS survey region. Unlike most 
spectroscopic surveys, GAMA samples the entire SDSS r- 
band selected galaxy population in an extensive and unbi- 
ased fashion to r pct = 19.8 mag - although for delivery of 
robust photo-z's to the full depth of SDSS, it is necessary to 
supplement GAMA with deeper data. Below we briefly sum- 
marize the method, explained in greater detail in Parkinson 
et al. (in prep). 

We use the Artificial Neural Network code ANNz 
(Collister & Lahav 2004) with a network architecture of 
N:2N:2N:1, where N=6 is the number of inputs to the net- 
work (5 photometric bands and one radius together with 
their errors). In order to develop a photo- z code valid for 
the full-SDSS region we use the ugriz iibercal calibration 
of SDSS DR7 (Padmanabhan al. 2008), together with the 
best of the De Vaucouleurs or Exponential half-light radius 
measurement for each object (Abazajian et al. 2009). Our 
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Figure 11. Original SDSS Pctrosian magnitudes versus our SEx- 
tractor based Kron magnitudes (top), original SDSS Petrosian 
magnitudes versus our GALFIT3 Sersic magnitudes (integrated 
to 10-Re) (middle), and our SExtractor Kron versus our GAL- 
FITS Sersic magnitudes (integrated to 10R e ) (lower). Contours 
vary from 10 per cent to 90 per cent in 10 per cent intervals. 



training set is composed of all GAMA spectroscopic redshifts 
and zCOSMOS spectroscopic redshift (Lilly et al. 2007), as 
well as the 30-band COSMOS photometric redshifts (Ilbert 
et al. 2009). With a redshift error of < 0.001, these are as 
good as spectroscopy when it comes to calibrating 5-band 
photo- z's from SDSS alone. The corresponding fractions in 
the training, validating and testing sets are 81 per cent, 4 
per cent and 15 per cent respectively, with over 120k galax- 
ies in total. For COSMOS (and hence zCOSMOS) galaxies, 
we perform a careful matching to SDSS DR7 objects. This 
large, complete, and representative photo-z training set al- 
lows us to use an empirical regression method to estimate 
the errors on individual photometric redshifts. The left panel 
of Fig. [l2]shows that up to z p h ~ 0.40 the bias in the median 
recovered photo-z is less than 0.005/(1 + 2 p h), while beyond 
it (not shown) the bias is typically within 0.01/(1 + 2 p h)- 
Considering the random errors on the individual estimates, 
this systematic bias is negligible, while still very well quan- 
tified. The right panel of Fig.[l2]presents the equivalent plot 
using SDSS DR7 photometric redshifts, as listed in the Pho- 
toz table of the SDSS Sky Server. In both cases, these plots 
show only galaxies with r pct < 19.8 mag. This comparison 
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Figure 13. Spectroscopic redshift distribution of GAMA galax- 
ies (black), the corresponding photo-z distribution of GAMA 
galaxies (red) and one thousand Monte-Carlo generated photo-2 
distributions, including the Gaussianly distributed photo-z error 
(green) . 



highlights the importance of using a complete and fully rep- 
resentative sample in the training of photometric redshifts. 

For an extensive and complete training set such as 
GAMA, it is possible to derive a simple and direct method 
for calculating photo-z errors. We estimate the photo-z ac- 
curacy in a rest-frame colour magnitude plane with several 
hundreds of objects per 0.1 xO.l mag bin down to r pc t = 19.8 
mag, and typically about fifty objects per 0.1 x 0.1 mag bin 
down to r pot ~ 21.5 mag. This method provides a robust 
normally distributed error for each individual object. We 
note that with this realistic photo-z error it is possible to cor- 
rectly recover the underlying spectroscopic redshift distribu- 
tion with our ANNz photometric redshifts as illustrated in 
Fig. 1131 It is interesting to note that, without this added er- 
ror, the raw photo-z histogram is actually narrower than the 
spectroscopic histogram. This makes sense, because photo-z 
calibration yields no bias in true z at given z p hot: it is thus 
inevitable that there is a bias in z p hot at given true z, which 
is strong wherever N(z) changes rapidly; this effect can be 
seen in operation in the ranges z < 0.05 and 0.3 < z < 0.5. 
GamaPhotoz contains both photometric redshifts and errors 
for all systems to r mo( jci < 21.5 mag and the code is cur- 
rently being applied across the full SDSS region. We plan to 
release our prescription for these improved all-SDSS photo- 
z's by the end of 2010. In the meantime, we are open to col- 
laborative use of the data: contact gama@gama-survey.org. 
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Figure 12. (left) GAMA photometric and spectroscopic redshift comparison as a function of the photometric redshift for GAMA 
galaxies, colour coded as function of galaxy (over-)density. The thick red central line shows the median, while the outer lines delineate 
the central 68 per cent of the distribution as function of z pfl . (right) Same as the left panel, but for the SDSS photometric redshift 
estimate, extracted from the Photoz table in SDSS DR7. 



4 BUILDING AND TESTING THE GAMA 
CORE CATALOGUES 

We can now combine the five catalogues described above 
to produce the GamaCore catalogue of ~ 1 million objects, 
which contains what we define as our core product to r < 22 
mag. Figs I14H16I uses this catalogue to illustrate the pro- 
gression of the redshift survey over the years through the 
build-up of the cone plots in the three distinct regions. As 
the survey has progressed the fidelity with which structure 
can be resolved and groups identified has increased signif- 
icantly. In the final maps the survey is over 98% complete 
with a near uniform spatial completeness for all three regions 
thanks to the "greedy" targeting algorithm (see Robotham 
et al., 2010). From GamaCore we now extract three science 
ready subsets: 

(1) GAMA data release 1 (GamaCoreDRl), which in- 
cludes the majority of spectroscopic data acquired in year 1 
with r pot < 19.0 mag and a central strip in G12 (S ±0.5°) to 
rpet < 19.8 mag and which forms our primary public data 
release product at this time. 

(2) GAMA Science (GamaCoreMainSurvey), which in- 
cludes all objects within our main selection limits and suit- 
able for scientific exploitation by the GAMA team and col- 
laborators. 

(3) GAMA-Atlas (GamaCoreAtlasSV), which includes 
all year 1 and year 2 data with reliable matches to the 
H- ATLAS science demonstration catalogue as described in 
Smith et al. (2010) and is currently being used for science 
exploitation by the H- ATLAS team. For further details and 
for the H- ATLAS data release see http : / /www . h-atlas . org 

In the following section we explore the properties of the 
two primary catalogues (GamaCoreDRl and GamaCoreMain- 
Survey) and do not discuss GamaAtlasSV any further. 



4.1 Survey completeness 

In Fig. [6] we showed the completeness of the GAMA red- 
shift survey only. This was given by the number of objects 
with secure redshifts divided by those targeted during the 



three year spectroscopic campaign. In Fig. [17] we show the 
equivalent completeness but now for the combined data as 
described in section [2781 for our two principle scientific cat- 
alogues (GamaCoreDRl and GamaCoreMainSurvey). This is 
given by the known redshifts in the GAMA regions divided 
by the known targets in the GAMA region. The blue regions 
show the pre-existing redshifts acquired before GAMA com- 
menced and is mainly dominated by the SDSS DR6 with a 
spectroscopic limit of r pe t < 17.77 mag (cf. Table[2}. Gama- 
CoreDrl is shown in green where the main focus in year 1 
was to survey objects brighter than r pe t = 19 mag with the 
exception of a central strip in the G12 region (motivated by 
the desire to assess the completeness function in preparation 
for year 2 and 3 observations). 

GamaCoreDRl is relatively complete to r pet < 19.0 mag 
and the residual spatial bias can be compensated for by us- 
ing the survey masks provided (see section ■ After year 
3 observations (GamaCoreMainSurvey) one can see 95 per 
cent or greater completeness in all apparent magnitude in- 
tervals with a smooth tapering off from 99 to 95 per cent 
completeness in the faintest 0.5mag interval. As all targets 
without a redshift have now been visually inspected this fall 
in completeness is real and obviously needs to be taken into 
consideration in any subsequent analysis. However GAMA 
remains the most complete survey published to date with 
the 2dFGRS achieving a mean survey completeness of 90 per 
cent to bj = 19.45 mag, SDSS 90 per cent to r pct = 17.77 
mag, and the MGC 96 per cent to B — 20 mag. 

4.2 Survey bias checks 

In Fig. [T7] we see a small but non-negligible completeness 
bias with apparent magnitude, hence it is worth checking 
for any bias in the obvious parameter space of colour, sur- 
face brightness, concentration, and close pairs. We explore 
the first three of these in Fig. [TS] for year 1 data (left), year 
2 data (centre) and year 3 data (right). The colour indi- 
cates the degree of completeness in these two-dimensional 
planes and the contours show the location of the bulk of the 
galaxy population. In all three cases there is no strong hori- 
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Figure 14. Redshift cone diagram for the GAMA 09 region showing (top to bottom), pre-existing data, year 1 data release added, year 
2 data release added, year 3 data release added. 



zontal bias that would indicate incompleteness with colour, 
surface brightness, or concentration. This is not particularly 
surprising given the high overall completeness of the sur- 
vey, which leaves little room for bias in the spectroscopic 
survey, but we acknowledge the caveat that bias in the in- 
put catalogue cannot be assessed without deeper imaging 
data. Note that in these plots surface brightness is defined 
by: /iefi = r pc t + 2.5 log 10 (2-7ra&) where a and b are the ma- 
jor and minor half-light radii as derived by GALFIT3 (see 
Kelvin et al. 2010 in prep, for details of the fitting process) 



and concentration is defined by log 10 (n) where n is the Sersic 
index derived in section [jT2l 

A major science priority of the GAMA programme is 
to explore close pairs and galaxy asymmetry (e.g., De Pro- 
pris et al. 2007). As GAMA is a multi-pass survey with each 
region of sky being included in 4-6 tiles (see. Fig. [2]), the 
close pair biases due to minimum fibre separations, which 
plague the 2dFGRS, SDSS, and MGC studies, are over- 
come. Fig. 1 191 shows the redshift completeness versus neigh- 
bour class. The neighbour class of any galaxy is defined (see 
Baldry et al. 2010) as the number of target galaxies within a 
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Figure 15. Redshift cone diagram for the GAMA 12 region showing (top to bottom), pre-existing data, year 1 data release added, year 
2 data release added, year 3 data release added. 



40" radius. The nominal 2dF fibre-collision radius. The fig- 
ures shows the progress towards resolving these complexes 
at the outset and after each year the survey has been op- 
erating. All low neighbour class objects are complete with 
only a few redshifts outstanding in a small fraction of the 
higher neighbour class complexes. The data therefore rep- 
resents an excellent starting point for determining merger 
rates via close pair analyses. 



4.3 Overdensity/underdensity of the GAMA 
regions 

All galaxy surveys inevitably suffer from cosmic variance, 
more correctly stated as sample variance. Even the final 
SDSS, the largest redshift survey to date, suffers an esti- 
mated residual 5 per cent cosmic variance to z < 0.1 (see 
Driver & Robotham 2010) — mainly attributable to the 
"Great Wall" (see Baugh et al. 2004; Nichol et al. 2006). 
As the GAMA regions lie fully within the SDSS it is possi- 
ble to determine whether the GAMA regions are overdense 
or underdense as compared to the full SDSS over the un- 



© 2010 RAS, MNRAS OOO.IWTl 



18 Driver et al. — DRAFT v3.1 



Pre -existing redshifts (GAM A 15) 

0.2 

0.1 

0.0 — 




Redshift 



Figure 16. Redshift cone diagram for the GAMA 15 region showing (top to bottom), pre-existing data, year 1 data release added, year 
2 data release added, year 3 data release added. 



biased redshift range in common (z < 0.1) and relative to 
each other at higher redshifts where the SDSS density drops 
due to incompleteness imposed by the r pe t < 17.77 mag 
spectroscopic limit. Fig. [501 shows results with the upper 
panel showing the density of M* ±1.0 galaxies (defined as 
M* — 5 log ft = -20.81 mag) in 0.01 redshift intervals, for 
a 5000 sq.deg. region of the SDSS DR7 (blue), and for the 
three GAMA blocks (as indicated) and the sum of the three 
regions (solid black histogram). The k and e corrections 
adopted and the methodology, including the SDSS sample 
construction, are described in Driver & Robotham (2010) 



but are not particularly important as we are not exploring 
trends with z but rather the variance at fixed z. The error 
bars shown are purely Poisson and so discrepancies larger 
than the error bars are indicative of cosmic variance induced 
by significant clustering along the line-of-sight. The middle 
and lower panels show the sum of these density fluctuations 
relative to either SDSS (middle) or the average of the three 
GAMA fields (lower). Note that these panels (middle, lower) 
now show the density fluctuation out to the specified redshift 
as opposed to at a specific redshift (top), i.e., the cosmic vari- 
ance. We can see that out to z = 0.1 the three GAMA fields 
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Figure 17. The final completeness of the combined redshift catalogues prior to survey commencing (blue), after GamaCoreDRl (green) 
and GamaCorcMainSurvey (red) for G09, G12 and G15 (left to right). In all cases the solid black lines shows 95 per cent completeness. 
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Figure 19. The redshift completeness as a function of Neighbour 
Class (NC). NC is defined for each target as the number of other 
Main Survey targets within 40" (numbers with each NC value 
are annotated in the figure). The various lines show the progress 
in resolving clustered objects at the outset and after each year of 
observations for the r-limited Main Survey. Note the bias toward 
NC>0 after Year 1 and 2 is the result of increasing the priority 
of clustered targets (Robotham et al. 2010) in order to avoid 
being biased against NOO after Year 3. The final survey has 
fully resolved almost all close complexes leaving only a minimal 
bias when determining merger rates via close pairs: the redshift 
completeness is > 95% for NC < 5. 



are overall 15 per cent under dense with respect to a 5000 
sq.deg. region of SDSS DR7 (NB: the volume surveyed by 
the SDSS comparison region is ~ 1.3 x 10 7 /i~ 3 Mpc 3 ). This 
is extremely close to the 15 per cent predicted from Table 2 
of Driver & Robotham (2010). Beyond z = 0.1 one can only 
compare internally between the three GAMA fields. The in- 
ference from Fig. [20] is that for any study at z < 0.2 the 
cosmic variance between the three regions is significant with 
G09 in particular being under-dense for all redshifts z < 0.2 
when compared to the other two regions (NB: the volume 
surveyed within the combined GAMA regions to z < 0.2 is 
~ 2.8 x 10 6 /i _3 Mpc 3 ). It is therefore important when consid- 
ering any density measurement from GAMA data to include 
the cosmic variance errors indicated in Fig. [20] For example 
a luminosity function measured from the G09 region only 
out to z — 0.1 would need to be scaled up by xl.41 etc. 



5 MASKS 

For accurate statistical analysis of GAMA it is essential 
to have a full understanding of the criteria that define 
its parent photometric catalogue, and also of the spatial 
and magnitude-dependent completeness of the redshift cat- 
alogue. For this purpose we have defined three functions 
characterizing this information as a function of position on 
the sky and magnitude selection. Here we present briefly the 
two most important ones, i.e. the survey imaging mask func- 
tion and the redshift completeness mask relative to the main 
r pe t survey limits. The combination of both functions is key 
for any spatially dependent measurements based on GAMA 
data, and in particular on GAMA year 1 data. Additional 
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Figure 18. Completeness is the bivariate planes of apparent magnitude and u — r colour (top), apparent magnitude and effective surface 
brightness (middle), apparent magnitude and Sersic index (bottom) for GAMA data after year 1 (left), year 2 (centre) and year 3 (right). 
Contours show the percentage of galaxies enclosed from out to in 99, 95, 75, 50 & 5 per cent. Apart from the progression in apparent 
magnitude over the three years no other obvious bias is evident. 



survey completeness masks for other selections (e.g., z, K) 
or any combination of are available to the Team and will be 
released shortly. 

5.1 The imaging mask 

We are interested in knowing which regions of the GAMA 
areas have not been properly covered by the SDSS imaging 
survey, or which should be excluded owing to the presence 
of bright stars. For that reason, we want to map out SDSS 
imaging areas containing any of the following information: 
bleeding pixel, bright star, satellite trail or hole. 

First we create, following the imaging mask information 
available on the SDSS DR6 websit«0, the associated convex 

3 http: //www. sdss . org/dr6/products/images/use_masks .html 
and http : //www. sdss . org/dr6/algorithms/masks .html 



polygons, delimiting areas for which imaging information is 
either not available or could be corrupted (as in the case of 
bleeding pixels). We primarily use the r-band imaging mask 
information, as GAMA is nominally an r-band selected sur- 
vey, but for completeness we construct all five SDSS imaging 
masks for the GAMA areas. 

Then we build an additional bright star mask based on 
stars down to V < 12 in the Tycho 2, Tycho 1 and Hipparcos 
catalogues. For each star we define an exclusion radius r, 
defined as: 

r = R s /0.8 for 10 < V < 12 (2) 

r = R s /0.5 for V < 10 (3) 

where R a is the scattered-light radius, as estimated based 
on the circular region over which the star flux per pixel is 
greater than five times the sky noise level. Further details 
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5.2 The r pct redshift completeness mask 
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Normally a simple way to define a redshift success rate would 
be to make use of the geometry defined by the complete set 
of 2° fields that were used to tile the survey region for spec- 
troscopic observations. However, due to the tiling strategy 
adopted in the first year, due to the preselection of observing 
the brightest targets first, and due to the much higher num- 
ber density of galaxies than can be accommodated in a sin- 
gle 2dF field, this simple and straightforward approach does 
not account well enough for the spatial incompleteness of the 
survey. Therefore for year 1 data, we had to develop different 
completeness masks each defined for a different magnitude 
limit interval: 



• For G09 and G15, there are two completeness masks 
each: one for galaxies brighter than r pct = 19.0 and one for 
19.0 < r-pot < 19.4. 

• For G12, there are three completeness masks: the same 
two as for G09 and G15, as well as one for 19.4 < r pot < 19.8. 



Figure 20. (upper) The differential number-density of M* ± 1.0 
mag galaxies in redshift intervals of 0.01. (middle) the overdensity 
of M* ±1.0 mag galaxies out to the specified redshift for the 
survey indicated by line colour relative to that seen in a 5000 
sq.deg. region of SDSS DR7. (lower) the overdensity of M t ± 1.0 
mag galaxies out to the specified redshift for the survey indicated 
relative to that seen over all three GAMA regions. 



on the bright star exclusion mask are given in section 3.3 of 
Baldry et al. (2010). 

The final imaging mask function is then the union of the 
two separate functions. For ease of use, the imaging masks 
have been pixelated using an equal area projection. 



Once the samples for which completeness masks are 
needed have been defined, one just needs to provide a rea- 
sonable definition for the redshift success rate. We choose to 
tessellate the GAMA regions with a large number of sectors. 
We incorporate the relevant imaging mask at this stage, by 
imposing that sectors do not cover any regions masked by 
the imaging mask. Each sector contains between 15 and 50 
galaxy targets, is limited in extent (less than 24 arcmin) 
and in size (less than 225 arcmin 2 ). These conditions are 
necessary to avoid shot-noise dominated masks, and to guar- 
antee that small scale information is preserved as much as 
possible. We note that in the current implementation these 
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sectors are not uniquely defined but once specified any 
given position on the sky belongs to a unique sector. For 
each sector, 6, we define the redshift success rate, R z {8), for 
a sample of galaxies within specified magnitude limits, as 
the ratio of the number of galaxies for which good quality 
redshifts have been obtained, N z {&), to the total number 
of objects contained in the tiling catalogue, Nt(6). The red- 
shift completeness of a given sector, R z (0), should be clearly 
distinguished from the redshift completeness of a given 2dF 
field, cf, since multiple overlapping fields can contribute to 
a single sector and cf is a measure of the quality of the 
observing conditions for galaxies observed at the same time. 

We present in Fig. the completeness masks for all 
three regions for GamaCoreDRl (left) and GamaCoreMain- 
Survey (right) to r pct < 19.0. Completeness masks for all 
regions and different selections are available to the team 
members and collaborators only at this stage. 



CO - 
C\J - 








c 
C 


"GAMA 
Non-GAMA 






DEC 

1 

,1,1, 





















t 1 1 1 1 1 ^ 1 r^— | 1 1 r 

130 132 134 136 138 140 
RA 



Figure 22. The location of the H-ATLAS SV region overlaid on 
the GAMA 9hr region (dotted lines) and showing the location 
of pre-existing (blue) and new GAMA (red) redshifts within the 
common region. 



6 GAMA SCIENCE READY CATALOGUES 
AND DATA RELEASE 1 

The combination of our five input catalogues, as outlined 
in section U constitutes our core GAMA catalogue of ~ 1 
million galaxies lying within the GAMA regions and extend- 
ing to approximately r pe t < 22.0 mag. This catalogue has 
inhomogeneous selection, is liable to be incomplete towards 
the faint-end, along with significant noise in the photometry 
at the very faint limit and spurious detections. We therefore 
extract from this dataset three science-ready catalogues and 
one overflow catalogue. These four catalogues along with the 
complete SWARP'ed mosaics of the GAMA regions, associ- 
ated spectra bundles, images, Sersic profiles and a variety 
of data inspection tools (including a MySQL access point) 
constitute our first data release and are now available via: 
: / / www . gama- survey. org/ 



6.1 GamaCoreMainSurvey - Available 01/07/12 

The GamaCoreMainSurvey is the GAMA team's principal sci- 
ence catalogue and is constructed from GamaCore by remov- 
ing all objects outside our Main Survey limits as defined in 
Baldry et al. (2010). These limits are: r pot < 19.4 mag in 
G09 & G15 and r pct < 19.8 mag in G12 (114 441 objects). 
Kktou < 17.6 with r modo i < 20.5 mag (61393 objects) or 
^Modei < 18-2 with r model < 20.5 mag (55534 objects) and 
are selected using SURVEY_CLASS> 3. For a description of 
the SURVEY_CLASS parameter see Table [3|. This amounts 
to 119 778 objects in total of which 101 576 are new redshifts 
provided by GAMA and 18202 pre-existing. Fig. [T7] shows 
the completeness in the r-band however similar plots can 
trivially be constructed in z or K. The spectroscopic com- 
pleteness (which includes objects not targeted) is 98.2 per 
cent in r, 99.3 per cent in z and 98.6 per cent in K. The pa- 
rameters contained in this catalogue are listed in Table lA"2l 
As this is our main science catalogue we place an embargo 
on its public release until l si July 2012 but are open to 
requests for collaboration sent to: gama@gama-survey.org 

4 This could be achieved by making each target galaxy the centre 
of a subsector and then increasing the radius of all subsectors and 
create sectors by the merger of overlapping subsectors. 



6.2 GamaCoreDRl — Available 2010-06-25 

The GamaCoreDRl is the current GAMA Data Release 1 
catalogue and is constructed from GamaCoreMainSurvey by 
applying a strict r-band selection of r pet < 19.4 mag in G09 
& G15 and r pe t < 19.8 mag in G12 (this is implemented by 
extracting objects with SURVEY.CLASS > 5. While the final 
GAMA survey is independently r, z and K selected, a strict 
r-band cut is imposed because observations during year 1 
were r-band limited only with the additional cuts coming 
in in years 2 and 3. Future releases involving year 2 and 
3 data will therefore include z and K band selected sam- 
ples. Selecting on SURVEYXLASS > 5 yields a catalogue 
of 114441 objects spread across the three GAMA regions. 
Table IA2I contains a list of the parameters being released 
at this time. Redshift information (i.e., quality, signal-to- 
noise etc.) are provided for all targets while the redshifts 
are only provided for year 1 observations above r pct = 19.0 
mag (52324 objects) and for a deeper narrow strip in G12 
with declination 5 = ±0.5° to 19.0 < r pot < 19.8 mag (7533 
objects). These limits are chosen as they represent our year 
1 spectroscopic targeting limits (see Figs 1171 1181 fc I2T1 left) . 
The redshift release constitutes approximately 50 per cent 
(59479 objects; 41902 GAMA; 17577 pre-existing), of the 
entire GAMA dataset (or 32.7 per cent of the redshifts ac- 
quired by the GAMA team over the past three years). This 
dataset provides a clean, well defined, r-band selected sam- 
ple suitable for scientific exploitation as shown by Fig. 1171 
and Fig. 1181 Please note that redshifts outside the speci- 
fied flux limits or collected after year 1 are currently de- 
noted as ZJ3ELIO = —2. The redshift quality and signal- 
to-noise information refers to the final best spectrum held 
in GamaCoreMainSurvey providing advance information as 
to whether a reliable redshift exists in the larger database. 
Requests for individual redshifts for specific objects can be 
directed to: gama@gama-survey.org 

6.3 GamaCoreAtlasSV — Available 2010-06-25 

The GamaCoreAtlasSV is the subset of objects that have 
both a GAMA and H-ATLAS detection (see Smith et al. 
2010 for details of source matching). 1175 objects coexist 
between the two data sets in a region defined by the geomet- 
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Table 3. The SURVEY_CLASS parameter given in GamaTiling. Objects take the higher value if satisfying more than one criteria. All objects 
with SURVEY_CLASS > 2 pass the standard star-galaxy separation. See Baldry et al. (2010) for details. 



value 


number 


frac. Q > 3 


criteria and notes 


7 


59756 


0.995 


r pet < 19.0 and A sg > 0.25 (high priority in Year 1) 


6 


54685 


0.968 


r'pet < 19.4 or r pe t < 19.8 in G12 (> 6 for r-limited Main Survey) 


5 


1386 


0.703 


^modcl < 18.2 (> 5 for SDSS-mag limited Main Survey) 


4 


3951 


0.879 


Jf AB < 17.6 (> 4 for Main Survey) 


3 


33134 


0.343 


19.4 < r pct < 19.8 in G09 or G15 (F2 fillers) 


2 


14861 


0.196 


Smodel < 20.6 or r modcl < 19.8 or i mo dcl < 19.4 in G12 (F3 fillers) 


1 


1255 


0.577 


radio selected targets (Fl fillers, Ching et al. in preparation) 





822 


0.394 


2 < VIS_CLASS < 4 or fails star-galaxy separation for r^j, < 17 



0.1 0.2 0.3 0.4 0.5 0.5 0.7 
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Figure 23. (Main panel). The absolute magnitude versus red- 
shift plane for the three science ready catalogues (as indicated), 
(left) the distribution collapsed in redshift, (upper) the distribu- 
tion collapsed in absolute magnitude. 

ric overlap between G09 and the H- ATLAS SV region (see 
Fig. l22[) . with an approximate area of ~ 16.7 deg 2 . Using the 
GamaCoreMainSurvey catalogue we find ~ 15 per cent of 
GAMA Main Survey targets have a H- ATLAS source match. 
Of these matches 89 do not have any available redshift, 306 
have SDSS redshifts, 792 have GAMA redshifts and 1 object 
has a 2SLAQ-QSO redshift. The GAMA redshifts comprise 
of 523 year 1 redshifts, 182 year 2 redshifts and 87 year 3 
redshifts. There are no magnitude cuts imposed for Gama- 
CoreAtlasSV and therefore 178 of these redshifts derive from 
GamaCoreExtra. Note that to extract a homogeneously se- 
lected sample one must trim this catalogue to r pet < 19.4 
mag. Table IA3I shows the parameters included in this cata- 
logue. 



6.4 GamaCoreExtraDRl — Available 2010-06-25 

The GamaCoreExtra catalogue contains information on all 
objects for which redshifts are known but lie outside 
the GamaCoreDRl selection limits. In effect GamaCoreEx- 
tra serves as a redshift dumping ground and when com- 



bined with GamaCoreMainSurvey or GamaCoreDRl provides 
a full record of all objects with known redshifts in the three 
GAMA regions. This will include pre-existing redshifts at 
faint magnitudes in the public domain (predominantly from 
QSO and LRG surveys), along with filler objects targeted 
during the GAMA campaign. The parameters included are 
identical to GamaCoreDRl and therefore also described in 
Table lA"2l GamaCoreExtraDRl is an identical copy of Gama- 
CoreExtra except with the GAMA redshifts embargoed un- 
til a later release. GamaCoreExtra is highly inhomogeneous 
and volatile as additional redshifts come to light and es- 
sentially forms the starting point for any extension to the 
GAMA survey. Currently it contains ~ 16k galaxies of which 
~ 0.5k are pre-existing and ~ 15.5k have been acquired dur- 
ing GAMA observations. Requests for access to individual 
redshifts or redshifts for specified sub-regions should be di- 
rected to: gama@gama-survey.org 



6.5 Future plans 

The three science ready catalogues described in the previous 
sections are shown in Fig. [23] in the Absolute r-band mag- 
nitude (Ai r — 5 log h) versus redshift (z) plane (main panel), 
along with the collapsed N(z) and N(M r ) distributions (up- 
per and right panels respectively). The absolute magnitudes 
are calculated using fc-corrections derived from the GAMA 
u to K photometry, and total throughput (including atmo- 
sphere) temperature adjusted filter transmission curves, and 
the KCORRECT software package (Blanton & Roweis 2007). 
The figure is intended to provide a rough indication of the 
parameter range being explored with L* galaxies routinely 
detected to z ~ 0.3 and objects as faint as M r — — 12 mag 
being recovered at z > 0.01. In future papers we will ex- 
plore more comprehensively the completeness of the initial 
input catalogue (GamaTiling) in the bivariate brightness dis- 
tribution plane (BBD; Hill et al. 2010, in prep.,) through 
a detailed BBD-analysis (see Driver et al. 2005), as well as 
through comparisons to deeper and higher spatial resolution 
data. 

In the near term the GAMA team will be produc- 
ing a number of extended data products or data manage- 
ment units (DMUs), which include: spectral line analysis via 
GANDALF/ppxf (Sarzi et al., 2006; Cappellari & Emsellem 
2004) and stellar population model fits (Maraston 2005; 
Thomas, Maraston & Bender 2003), bulge-disc decompo- 
sition via GALFIT3, stellar masses and revised photometric- 
redshifts via INTEREST, environmental markers, and group 
catalogues. These data products will be outlined in more 
detail in our next data release later in the year (Liske et 
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al., 2010, in prep.). Within the next 3 — 5 years we expect 
to ingest additional datasets from GALEX, HERSCHEL, 
WISE, GMRT, VISTA, VST, and ASKAP as well as ex- 
pand our core AAT programme to cover a larger area of 
~ 400 sq deg. to rsersic < 19-8 mag. In due course DMUs 
on multiwavelength SED and HI modeling will provide dust 
mass/temperature measurements, and gas and dynamical 
mass measurements. 



7 SUMMARY 

The GAMA basic-data including GAMA and pre-GAMA 
spectra and redshifts, GAMA SWARP processed images 
in ugrizYJHK, u to K matched aperture photometry, 
corrections to total magnitudes are now available at: 
http://www.gama-survey.org/. GAMA Data Release 1 
(GamaCoreDRl) data is publicly available via a downloads 
page immediately, and the full catalogue (GamaCoreMain- 
Survey) is available for use by the GAMA team and collab- 
orators. In addition we make publicly available a catalogue 
(GamaCoreAtlasSV) including year-2 data for the H-Atlas 
Science Verification region, which form the basis for a num- 
ber of H-Atlas Science Verification papers. 

Data inspections tools, including a MySQL tool and di- 
rect data inspection toolkit, are provided and will be devel- 
oped further in due course. The sum of year 1, 2, and 3 data 
constitutes GAMA Phase I, which is now complete having 
utilised archival data from the Sloan Digital Sky Survey, the 
UK Infrared Deep Sky Survey, and the NASA Extragalactic 
Database and an initial allocation of 66 nights of Anglo- 
Australian Telescope time. Value added data products will 
be released from January 1st 2011 and regularly thereafter. 
GAMA Phase II is now underway and aims to add in ad- 
ditional area coverage as well as to extend the depth of the 
survey. The community is encouraged and invited to con- 
tact the team for early access to GamaCoreMainSurvey. The 
contents of this catalogue are described in Tabic IA2I 



Group, Johns Hopkins University, the Joint Institute for 
Nuclear Astrophysics, the Kavli Institute for Particle As- 
trophysics and Cosmology, the Korean Scientist Group, the 
Chinese Academy of Sciences (LAMOST), Los Alamos Na- 
tional Laboratory, the Max-Planck-Institute for Astronomy 
(MPIA), the Max-Planck-Institute for Astrophysics (MPA), 
New Mexico State University, Ohio State University, Uni- 
versity of Pittsburgh, University of Portsmouth, Princeton 
University, the United States Naval Observatory, and the 
University of Washington. 

The UKIDSS project is defined in Lawrence et al. 
(2007). UKIDSS uses the UKIRT Wide Field Camera (WF- 
CAM; Casali et al, 2007). The photometric system is de- 
scribed in Hewett et al. (2006), and the calibration is de- 
scribed in Hodgkin et al. (2009). The pipeline processing 
and science archive are described in Irwin et al. (2009, in 
prep) and Hambly et al. (2008). We have used data from 
the 4th data release. 

This research has made use of the NASA/IPAC Ex- 
tragalactic Database (NED), which is operated by the Jet 
Propulsion Laboratory, California Institute of Technology, 
under contract with the National Aeronautics and Space 
Administration. 

The Millennium Galaxy Catalogue consists of imag- 
ing data from the Isaac Newton Telescope and spectro- 
scopic data from the Anglo Australian Telescope, the 
ANU 2.3m, the ESO New Technology Telescope, the Tele- 
scopio Nazionale Galileo and the Gemini North Tele- 
scope. The survey has been supported through grants 
from the Particle Physics and Astronomy Research Coun- 
cil (UK) and the Australian Research Council (AUS). 
The data and data products are publicly available from 
http://www.eso.org/~jliske/mgc/ or on request from J. 
Liske or S.P. Driver. 

Finally we would like to acknowledge the use of Topcat, 
Stilts, SExtractor, SWARP, PSFex, GALFIT3 and ALADIN 
astronomical software packages and the r data analysis pack- 
age. 
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APPENDIX A: GAMA CATALOGUES 

The four principle GAMA catalogues formed at this time 
are: GamaTiling, GamaCoreMainSurvey, GamaCoreDRl and 
GamaCore AtlasSv. The parameters contained in these cata- 
logues and their explanations are showing in Tables. IAD — 1A3I 
respectively. 
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Table Al. Parameters held in GamaTiling 



Column 


Parameter 


Units 


Definition 


1 


GAMAJD 


N/A 


Unique six digit GAMA identifier linked to GAMAs original SDSS DR6 input cat 


2 


SDSS.ID 


N/A 


Sloan Identifier from SDSS DR6 


3 


RA.J2000 


Degrees 


RA taken from SDSS DR6 


■1 


DEC.J2000 


Degrees 


Declination taken from SDSS DR6 


5 


r_FIBREMAG 


AB mag 


Fibre magnitude 


6 


r_PETRO 


AB mag 


Extinction corrected r-band Petrosian flux derived from SDSS DR6 


7 


uJVIODEL 


AB mag 


Extinction corrected model u mag from SDSS DR6 


8 


g_MODEL 


AB mag 


Extinction corrected model g mag from SDSS DR6 


9 


rJVIODEL 


AB mag 


Extinction corrected model r mag from SDSS DR6 


10 


LMODEL 


AB mag 


Extinction corrected model i mag from SDSS DR6 


11 


zJVIODEL 


AB mag 


Extinction corrected model z mag from SDSS DR6 


12 


NUM.GAMA.SPEC 


N/A 


Number of spectra available for this object 


13 


r.SB 


AB mag arcsec -2 


SDSS Pctro half-light r-band surface brightness 


11 


SG.SEP 


mag 


SDSS r_psf-r_model star-galaxy separation parameter 


15 


SG.SEP.JK 


mag 


GAMA AUTO J-K - f Jocus star-galaxy separation parameter 


16 


K.KRON.SELECTION 


AB mag 


Extinction corrected K magnitude used for selection (Old) 


17 


TARGET.FLAGS 


bit 


bitwise target criteria (see gama website) 


18 


SURVEY.CLASS 


N/A 


Survey class (see Table [3} 


19 


PRIORITY.CLASS 


N/A 


Priority class (see gama website) 


20 


NEIGHBOUR-CLASS 


N/A 


Number of Main Survey neighbours within 40 arcseconds (see Fig. 11 Q It 


21 


MASKJC.10 


N/A 


Mask value 0.0 to 1.0 around V < 10 mag stars (see Baldry et al. 2010) 


22 


MASKJC.12 


N/A 


Mask value 0.0 to 1.0 around V < 12 mag stars (see Baldry et al. 2010) 


23 


VIS.CLASS 


N/A 


Visual classification (sec Baldry et al. 2010) 


24 


VIS.CLASS.USER 


N/A 


Initials of visual classifier for VIS_CLASS 
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Table A2. Parameters held in GamaCoreMainSurvey, GamaCoreDRl, GamaCoreExtraDRl. 



lumn 


Parameter 


Units 


Definition 


1 


GAMAJAUJD 


N/A 


Unique IAU identifier 


2 


GAMAJD 


N/A 


Unique six digit GAMA identifier linked to GAMAs original SDSS DR6 input cat 


3 


SDSS.ID 


N/A 


Sloan Identifier from SDSS DR6 


■1 


RA.J2000 


Dej 


,rccs 


RA taken from SDSS DR6 


5 


DEC.J2000 


Dej 


;rees 


Declination taken from SDSS DR6 


6 


r_PETRO 


AB 


mag 


Extinction corrected r-band Pctrosian flux derived from SDSS DR6 


7 


Z.HELIO 


N/A 


Heliocentric rcdshift (—2 = embargoed. = —0.9 if nQ=l) 


8 


Z.QUALITY 


N/A 


Confidence on rcdshift (Definite, nQ=4, Reliable nQ=3, Uncertain nQ=2, Unknown 










nQ=l, Not observed nQ=99) 


9 


Z.SOURCE 


N/A 


Origin of redshift 1=SDSS DR6, 2=2dFGRS, 3=MGC, 4=2SLAQ-LRG, 5=GAMA, 










6=6dFGS, 7=UZC, 8=2QZ, 9=2SLAQ-QSO, 10=NED 


10 


Z.DATE 


N/A 


Date of observation if z_SOURCE=5 (i.e., GAMA), otherwise 99 


11 


Z.SN 


N/A 


Mean s/n of spectrum if z_SOURCE=5 (i.e., GAMA), otherwise 99 


12 


ZJD 


N/A 


Filename of best available spectrum. 


13 


PHOT.SOURCE 


N/A 


Photometric source (rd = r-band defined, sd = self-defined, sec Hill ct al. 2010a) 


14 


U.KRON 


AB 


mag 


Extinction u-band Kron magnitude derived from r-band aperture 


15 


g-KRON 


AB 


mag 


Extinction g-band Kron magnitude derived from r-band aperture 


16 


r_KRON 


AB 


mag 


Extinction r-band Kron magnitude derived from r-band aperture 


17 


i-KRON 


AB 


mag 


Extinction j-band Kron magnitude derived from r-band aperture 


18 


z.KRON 


AB 


mag 


Extinction z-band Kron magnitude derived from r-band aperture 


19 


Y.KRON 


AB 


mag 


Extinction Y-band Kron magnitude derived from r-band aperture 


20 


J.KRON 


AB 


mag 


Extinction J-band Kron magnitude derived from r-band aperture 


21 


H.KRON 


AB 


mag 


Extinction H-band Kron magnitude derived from r-band aperture 


22 


K.KRON 


AB 


mag 


Extinction X-band Kron magnitude derived from r-band aperture 


23 


U.KRON.ERR 


AB 


mag 


M-band Kron magnitude error 


21 


g_KRON_ERR 


AB 


mag 


(/-band Kron magnitude error 


25 


r_KRON_ERR 


AB 


mag 


r-band Kron magnitude error 


26 


LKRON.ERR 


AB 


mag 


•i-band Kron magnitude error 


2V 


Z.KRON.ERR 


AB 


mag 


z-band Kron magnitude error 


28 


Y_KRON_ERR 


AB 


mag 


Y-band Kron magnitude error 


29 


j.kron_err 


AB 


mag 


J-band Kron magnitude error 


30 


H_KRON_ERR 


AB 


mag 


H-band Kron magnitude error 


31 


K.KRON.ERR 


AB 


mag 


X-band Kron magnitude error 


32 


EXTINCTION^ 


AB 


mag 


Galactic magnitude extinction in r-band 


33 


r_SERS.MAG.10RE 


AB 


mag 


r-band Sersic magnitude truncated at 10 half light radii 
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Table A3. Parameters held in GamaCoreAtlasSV 



Column 


Parameter 


Units 


Definition 


1 


GAMAJAUJD 


N/A 


IAU certified GAMA ID 


2 


GAMAJD 


N/A 


Unique six digit GAMA identifier linked to GAMAs original SDSS DR6 input cat 


3 


SDSSJD 


N/A 


Sloan Identifier from SDSS DR6 


■1 


RA.J2000 


Dej 


;rccs 


RA taken from SDSS DR6 


5 


DEC.J2000 


Dej 


;rees 


Declination taken from SDSS DR6 


6 


r_PETRO 


AB 


mag 


Extinction corrected r-band Petrosian flux derived from SDSS DR6 


7 


Z.HELIO 


N/A 


Heliocentric redshift (—0.9 if nQ=l) 


8 


Z.QUALITY 


N/A 


Confidence on redshift (Definite, nQ=4, Reliable nQ=3, Uncertain nQ=2, Unknown 










nQ=l, Not observed nQ=99) 


9 


Z.SOURCE 


N/A 


Origin of redshift 1=SDSS DR6, 2=2dFGRS, 3=MGC, 4=2SLAQ-LRG, 5=GAMA, 










6=6dFGS, 7=UZC, 8=2QZ, 9=2SLAQ-QSO, 10=NED 


10 


Z.DATE 


N/A 


Date of observation if z_SOURCE=5 (i.e., GAMA), otherwise 99 


11 


Z.SN 


N/A 


Mean s/n of spectrum if z_SOURCE=5 (i.e., GAMA), otherwise 99 


12 


ZJD 


N/A 


Filename of best available spectrum. 


13 


PHOT.SOURCE 


N/A 


Photometric source (rd = r-band defined, sd = self-defined, see Hill ct al. 2010a) 


11 


U.KRON 


AB 


mag 


Extinction u-band Kron magnitude derived from r-band aperture 


15 


g_KRON 


AB 


mag 


Extinction g-band Kron magnitude derived from r-band aperture 


16 


r_KRON 


AB 


mag 


Extinction r-band Kron magnitude derived from r-band aperture 


17 


i-KRON 


AB 


mag 


Extinction j-band Kron magnitude derived from r-band aperture 


18 


z.KRON 


AB 


mag 


Extinction z-band Kron magnitude derived from r-band aperture 


19 


Y.KRON 


AB 


mag 


Extinction Y-band Kron magnitude derived from r-band aperture 


20 


J.KRON 


AB 


mag 


Extinction J-band Kron magnitude derived from r-band aperture 


21 


H.KRON 


AB 


mag 


Extinction //-band Kron magnitude derived from r-band aperture 


22 


K.KRON 


AB 


mag 


Extinction A'-band Kron magnitude derived from r-band aperture 


23 


U.KRON.ERR 


AB 


mag 


u-band Kron magnitude error 


24 


g_KRON_ERR 


AB 


mag 


(/-band Kron magnitude error 


25 


r_KRON_ERR 


AB 


mag 


r-band Kron magnitude error 


26 


i_KRON_ERR 


AB 


mag 


i-band Kron magnitude error 


2V 


Z.KRON.ERR 


AB 


mag 


z-band Kron magnitude error 


28 


Y_KRON_ERR 


AB 


mag 


Y-band Kron magnitude error 


29 


J-KRON.ERR 


AB 


mag 


J-band Kron magnitude error 


30 


h_kron_err 


AB 


mag 


//-band Kron magnitude error 


31 


K.KRON.ERR 


AB 


mag 


X-band Kron magnitude error 


32 


EXTINCTIONS 


AB 


mag 


Galactic magnitude extinction in r-band 


33 


r_SERS.MAG.10RE 


AB 


mag 


r-band Sersic magnitude truncated at 10 half light radii 


34 


HATLAS.IAU.ID 


N/A 


H-ATLAS ID as specified by Smith et al. (2010) 
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